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ABSTRACT 

We present a multi- wavelength study of the local Leo cold cloud (LLCC), a very nearby, very cold 
cloud in the interstellar medium. Through stellar absorption studies we find that the LLCC is between 
11.3 pc and 24.3 pc away, making it the closest known cold neutral medium cloud and well within the 
boundaries of the local cavity. Observations of the cloud in the 21-cm Hi line reveal that the LLCC 
is very cold, with temperatures ranging from 15 K to 30 K, and is best fit with a model composed 
of two colliding components. The cloud has associated 100 micron thermal dust emission, pointing 
to a somewhat low dust-to-gas ratio of 48 xlO -22 MJy sr _1 cm 2 . We find that the LLCC is too 
far away to be generated by the collision among the nearby complex of local interstellar clouds, but 
that the small relative velocities indicate that the LLCC is somehow related to these clouds. We 
use the LLCC to conduct a shadowing experiment in 1/4 keV X-rays, allowing us to differentiate 
between different possible origins for the observed soft X-ray background. We find that a local hot 
bubble model alone cannot account for the low-latitude soft X-ray background, but that isotropic 
emission from solar wind charge exchange does reproduce our data. In a combined local hot bubble 
and solar wind charge exchange scenario, we rule out emission from a local hot bubble with an 1/4 
keV emissivity greater than 1.1 Snowdens / pc at 3 a, 4 times lower than previous estimates. This 
result dramatically changes our perspective on our local interstellar medium. 
Subject headings: Galaxy: solar neighborhood ; ISM: HI; ISM: dust, extinction ; ISM: clouds 



1. INTRODUCTION 

The Galactic interstellar medium (ISM) is the perva- 
sive gas that fills our Galaxy. Our ISM has a huge range 
of temperatures and densities, from molecular gas at den- 
sities above 10 4 cm -3 and temperatures of only 10 K, to 
the hot ionized medium (HIM), which can exceed 10 6 
K in temperature, with densities dropping below 10~ 2 
cm~ 2 . The dynamic interplay between these different 
phases of the interstellar medium is not very well under- 
stood: How is the cold neutral medium (CNM) formed? 
What kind of gas fills low-density volumes in the ISM? 
How is it created? One way to approach these questions 
is to examine the ISM nearest the sun, an area called 
the local cavity, where our observations can most easily 
disentangle individual clouds, and examine their forma- 
tion, structure, interaction, and destruction. Indeed, the 
canonical basis for our understandi ng of the interplay 
betwe en hot and cold ISM phases, Mc Kee fc Ostriker] 
( 1977), takes observations of the local cavity as its start- 



ing point. 

Up until recently, the only ISM clouds thought to ex- 
ist near the sun were the complex of local interstellar 
clouds (CLIC), a group of about 15 low column den- 
sity (Nh ~ 10 18 cm -2 ), partial ly ionized clouds within 
about 15 pc of the sun (e. g. jRcdficld & Linsky |2008| 
RL08). Recent observations by Meyer et al.| (|2006|) have 



1 Department of Astronomy, Columbia University, Pupin 
Physics Laboratories, 550 West 120th Street, New York, New 
York, 10027. goldston@gmail.com 

2 Department of Astronomy, University of California, Berke- 
ley Berkeley, CA 94720 



Center for Interdisciplinary Exploration and Research in As- 



trophysics, Department of Physics and Astronomy, Northwestern 
University, Evanston, IL 60208 

4 Department of Physics and Astronomy, University of 
Louisville, Louisville, KY 40292 



added an entirely new animal to the local menagerie: a 
very cold CNM cloud, which we call the Local Leo Cold 
Cloud (LLCC). 

The first obser vation of the LLCC was conducted by 
Verschuur ( 1969[ ), who discovered a very cold component 



ot the ISM while looking towards intermediate velocity 
clouds. The very narrow observed linewidth in the 21-cm 
hyperfine transition of hydrogen led to his determination 
that this p arcel of ISM must have a t emperature of less 
than 30 K. Knapp & Verschuur ( 1972 ) followed this work 
up by mapping the LLCC (roughly a = 10 h, S = 10°), 
showing it to be composed of two angularly distinct com- 
ponents, eac h a few square degrees in area. |Crovisier fc| 



Kazes ( 1980[ ) found a spin temperature of about 20 K us- 
ing the Hi absorption-line towards the quasars 3C 225a 
and 3C 225b, but found no evidence for molecules in the 
1665 and 1667 MHz line of hydroxyl (OH). The LLCC 
was la ter independently discovered by |Heiles fc Troland| 
(2003), who saw it again in absorption toward 3 sources 
in their Arecibo millenniu m survey, and found tempera- 
tures of 14, 17 and 22 K. |Meyer et aL] ( |2006[ ) observed 
stars toward this set of clouds and found that every star 
in the direction of the cloud showed significant Na I ab- 
sorption at exactly the cloud's velocity in Hi, putting an 
upper limit of the dist ance to the cloud of 42 pc (HD 
83683). Most recently [Haud] pfllu] HaudlO) found the 
LLCC in the Leiden- Argentina-Bonn survey data, and 
found that it was connected to a much larger ribbon of 
clouds with consistent linewidths and a continuous veloc- 
ity distribution, stretching from -10° to +40° declination 
and from 7h to llh right ascension. This local ribbon of 
cold clouds (LRCC) extends through the constellations 
Sextans, Leo, Cancer, and Lynx, with small pieces in 
Hydra, Gemini, and Auriga. 
The discovery that there is an extended group of 
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nearby CNM clouds is intriguing in the context of the 
local cavity. The sun lives roughly in the middle of a 
largely evacuated volume of gas, roughly 200 pc wide in 
the Galactic plane, and perhaps somewhat pr olate (or 



even open) towards high Gala ctic latitude (e.g., Vergely 
eTaLp0l0| | Welsh et al |20lo| . The edges of this volume 



of gas are typically defined by the distance at which the 
total hydrogen column density reaches 10 19 3 cm -2 . The 
cavity seems to have a rather discrete edge, with very lit- 
tle material inside it. Since t he original Wisconsin rocket 
sounding experiments (e.g., IBurstein et al. 19771 it had 
been proposed that this cavity is tilled with an approxi- 
mately million-degree plasma, emitting thermally in soft 
X-rays. The million-degree plasma would explain the ob- 
served soft X-rays, which would be absorbed by neutral 
gas outside the cavity if its source were further away. 
Morphologically this made good sense as well; the bub- 
ble could have been carved out by a series of supernova 
explosions and stellar winds, wh ich would have supplied 
the energy to heat the gas (e.g., Cox & Reynolds 1987). 
This is in excellent agreement with theories of the ISM 
as a whole, with supernovae peppering the disk, blowing 
bubb les and launching material above the Galactic plane 
(e.g. , |de Avillez||2000| ). Once the ROSA T all-sky survey 



soft X -ray background ( RASS; SXRB; Snowdcn ^et al. 



19971 maps were available, Kuntz & Snowden ( 20001 pro 
posed a modified mechanism to explain the soft X-ray 
signature; part of the diffuse soft X-ray emission came 
from the halo and was partially absorbed by low col- 
umn, high-latitude Hi, and part of it originated from the 
local hot bubble. 

This picture has been challenged by a number of dif- 
ferent observations and theories, including concerns re- 
garding O vi interstellar absorption and the rel ationship 
between cavit y depth and soft X-ray flux (see |Welsh fe~| 
Shelton 2009 for a detailed summary), but none so con- 
cerning as the solar wind charge exchange (SWCX), a 
process that may be an alternate origin for the SXRB. 
ROSAT observations of comet Hyakutake showed signif- 
icant, unexplained X-rays, which were later determined 
to be from the electronic cascade of exchanged electrons 
from the relativ ely neutral come t to highly i onized so- 



lar wind metals ( Lisse et al.|1996 ). Both Cox (19981 and 
Freyberg ( 1998 ) suggested that this same SWCX process 
could generate X-rays in the heliopause, where the ion- 
ized solar wind met the largely neutral ISM of the local 
interstellar cloud (LIC). Accurate models of the SWCX 
emission are very difficult to generate, as the efficiency 
of exchange is largely unknown and the atomic physics 
of high ions is not fully constrained; as such there is no 
clear consensus as to how much of the soft X-ray emis- 
sion could be coming from SWCX, rather than a LHB . 
Tha t being said, models by|Robertson fe Cravens| ( 2003 ) 
and Koutroumpa et al. ( 200U ) show that somewhere be- 
tween 50% and 100% of the X-rays could be originating 
from SWCX, which leaves the fate of the LHB X-rays 
very uncertain. 

It has been suggested by various auth ors (e.g., Cox & 
Reynolds||T987l |Welsh fe Shelton| |2009| that an excel- 
lent way to determine the provenance of the soft X-rays 
would be to find a nearby cloud capable of providing a 
screen to soft X-rays. If a shadow for this cloud could 
be found, we could distinguish between X-rays coming 
from a LHB and X-rays coming from SWCX. We believe 



that the LLCC is just such a cloud, at the right distance 
and column density, and we devote a large fraction of 
this paper to carefully determining the crucial parame- 
ters needed to use the LLCC in this proposed shadowing 
experiment. 

In Section [2] we describe the observation and data re- 
duction of a multi- wavelength data set, including new 
observations in the optical and radio, as well as archival 
data in the infrared and X-ray. In Section [3] we discuss 
how we analyze these data in the context of the LLCC. 
In Section [4] we discuss various properties of the cloud 
we can determine from the observations. In Section [5] 
we discuss the implications of the LLCC observations for 
the contents of the local cavity, including the relation- 
ship to CLIC clouds and the implications for the LHB 
and SWCX theories. We conclude in Section [6j 

2. OBSERVATIONS AND DATA REDUCTION 

One of the main advantages of studying such a nearby 
structure is its scale on the sky; HaudlO detects the 
LRCC over 86 square degrees, and we find the LLCC 
alone to cover 22 square degrees. Because of this large 
area, we are capable of detailed morphological studies 
with low resolution observations. In particular, we can 
use archival all-sky surveys to pin down the structure 
of the object over a broad range of wavelengths. Ad- 
ditionally, the large area allows us to study the object 
more easily in absorption as more background sources 
are available. 

2.1. Radio: GALFA-HI 

The LLCC was originally discovered in the 21-cm (1420 
MHz) hyperfine transition of Hi, and this kind of obser- 
vation remains the best way to study the morphology of 
the structure. Temperat ure measurements of ^20 K fr om 
Crovisier fe Kazes| ( |1980] ) and |Heiles fe Troland| ( |2003| ) in- 
dicate that the cloud should be largely neutral. Measured 
column densities below 10 21 cm -2 and the non-detection 



of OH (|Crovisier fc Kaz cs 1980) indicate that the cloud 
should be largely atomic and not entirely optically thick 
to Hi. Thus, Hi is an excellent tracer of the bulk of the 
cloud. 

We observed the LLCC as part of the Galactic Arecibo 
L-Band Feed Array Hi (GALFA-Hi) survey, an ongo- 
ing survey of Hi conducted with the Arecibo 305m tele- 
scope using the Arecibo L-band Feed Array (ALFA). The 
GALFA-Hi survey is a survey of all Galactic Hi within 
the Arecibo observing range (—1° < S < 38°), at a spec- 
tral resolution of 0.184 km s _1 , and an angular resolution 
near 4'. The details of the observation methods, data re- 
duction, and properties of the GALFA-Hi survey can be 
found in Peek et al. (2011). Observations for this region 
were conducted both as part of a targeted map, and as 
part of the Turn On GALFA Survey, a commensal (si- 
multaneous) survey conducted alongside the the Arecibo 
legacy fast ALFA project. The targeted map was con- 
ducted in a high-speed, basketweave or meridian-nodding 
mode, whereas the Turn On GALFA Survey data were 
collected in a low-speed drift mode. Most areas of the 
map were covered by both modes, and have an RMS 
noise of 0.17 K over a 0.184 km s _1 channel, but a few 
areas were only covered with the basketweave mode, and 
have an RMS noise of 0.36 K per 0.184 km s _1 channel. 
It is important to point out that the detailed analysis of 
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the Hi column density found in §3. 1| could not be accom- 
plished without the exquisite spectral resolution of the 
GALFA-Hi survey. We note that while we refer to the 
section of cloud we are examining as the Local Leo Cold 
Cloud, the southern component of the LLCC is actually 
in Sextans. 

2.2. Infrared: IRAS & IRIS 

A key result from the Infrared Astronomical Satellite 
(IRAS) was the discovery of the correlation between the 
infrared diffuse emission and the Galac tic Hi column 

[1984} 
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density for |6| > 10° ( Low et al 
sion reached by these authors, and 
is a relatively uniform interstellar radiation field in the 
Galaxy and a relatively constant large-grain dust frac- 
tion in the diffuse Hi. The interstellar radiation field 
heats the dust, which reradiates in the infrared, which 
in turn creates the observed correlation between Hi col- 
umn density and dust emission. Thus, the investigation 
of the large-grain dust contents of the LLCC can be con- 
ducted by examining the infrared flux in the IRAS data. 
IRAS observed the vast bulk of the sky for 300 days 
in 1983, in the infrared wavelength band s centered on 
12, 25, 60, and 100 ^m ( |Beichman||l987} . These data 
have been reduced in numerous ways and have produced 
many different maps: the original SkyFlux Atlas, the 



IRAS Sky S urvey Atlas, the maps produced in [Schlcgel 
" ille 



et al. ( 1998), and the IRIS maps described in Mivilie 
Deschenes & Lagache ( 2006[ ). We use these most recent 
IRIS maps, as they are produced at 4' resolution, nearly 
identical to the resolution of the GALFA-Hi maps. In ad- 
dition to the IRIS reduction, we further clean the data 
set by removing point sources down to a magnitude 1 Jy. 
We remove these point sources by fitting with Gaussian 
PSF, and excise any point sources we cannot cleanly re- 
move. We note that the presence of IRC+10216 in this 
field has made the analysis of the 60, 25, and 12 mi- 
cron band impossible, as this exceedingly bright star has 
contaminated a large swath of the map. 

2.3. Optical: KPNO & CPS 

Based on observations of the interstellar Nai D2 
A5889.951 and Di A5895. 924 absorption towa rd 33 stars 
in the LLCC sky region, Meyer et al. (2006) were able 
to place a significant upper limit on the L LCC distance. 
Using the updated Hipparcos parallaxes of |van Leeuwen| 
(2007), this upper limit corresponds to the distances 
(39.91^3 and 40.5±J;J pc) of the nearest stars (HD 83808 
and HD 83683) exhibiting LLCC Nai absorption. The 
challenge in tightening this distance constraint further is 
finding nearer stars within the LLCC Hi 21-cm emission 
sky envelope that are bright enough for high-resolution 
optical absorption-line spectroscopy. Given the paucity 
of such stars, another possibility is to consider nearer 
stars just outside the Hi envelope that might be bright 
enough to detect very weak optical absorption at a gas 
column threshold more sensitive than the 21-cm obser- 
vations. 

Fortunately, there is a very bright (V = 1.35) star, 



van 



Regulus (a Le o), at a distance of 24.3 ± 0.2 pc 
Leeuwen||2007 ) and a s ky position that is 3° outside the 
LLCC 21-cm envelope. Meyer et al. ( 2006 ) obtained high 
S/N spectra of Regulus in the JN a 1 wavelength region to 



divide out the atmospheric absorption lines in the spec- 
tra of their program stars. These Regulus spectra reveal 
no evidence of any interstellar Nai absorption to high 
precision. Since Regulus is often used as a spectral stan- 
dard for observations in other wavelength regions, we 
examined our archive of data obtained with the 0.9 m 
coude feed telescope and spectrograph at Kitt Peak Na- 
tional Observatory (KPNO) over the past fifteen years. 
Observations of Regulus in the vicinity of the interstel- 
lar Ca II K A3933.663 and H A3968.468 lines were ob- 
tained with this instrumentation in November 2002 in 
support of a study of variable interstella r absorption to- 
ward p Leo ( |Lauroesch fc Meyer|[2003| . Eight F3KB 
CCD echelle spectra spanning a total exposure time of 
80 minutes were taken of Regulus with a spectrograph 
configuration yielding a measured velocity resolution of 
3.6 km s _1 at the Ca II K wavelength. We utilized the 
NOAO IRAF echelle data reduction package to bias- 
correct, scattered light-correct, flat-field, wavelength- 
calibrate, order-extract, and sum the individual CCD 
exposures into the final net Regulus spectrum. 

For the lower limit on the LLCC distance, we bor- 
rowed archival data from the California Planet Search 
of the nearby M-dwarf HI P 47513. This star is at a dis- 
tance of 11.26 ± 0.21 pc ( |van Leeuwen 2007). In addi- 
tion to HIP 47513, we used spectra of tour additional M- 
dwarfs with similar features as comparison stars. Those 
other spectra allow us to compare the Nai Di line core 
of HIP 47513 with the cores of other M-dwarfs, and thus 
to place a stronger constraint on the interstellar Na 1 col- 
umn along the HIP 47513 line of sight. 

Our five M-dwarf spectra were collected with HIRES, 
the high-re solution echelle spectrograph on Keck I ( Vogt 
et al.|1994 1 , between 2005 February and 2006 December. 
Details are in Table [I] The spectra were originally col- 
lec ted with the intention of detecting extrasolar planets; 
seel Marcy fc Butler|1992] for a full description of the CPS 
and its goals. Because of their radial velocity measure- 
ment technique, the vast CPS archive contained only a 
few useful spectra of each star. The CPS places an iodine 
cell in the light path of its stars, which creates a forest 
of I2 absorption lines for radial velocity reference. Here 
we make use only of the iodine-free "template" spectra, 
which contain no I2. 

The spectra have a resolution R — 70, 000 at A = 
5500A. We chose the stars in Table [l] because they were 
M-dwarfs with B - V within ±0.12 of HIP 47513, and 
they had iodine-free observations in the CPS archive. 
We began with 33 candidate stars that fit those two cri- 
teria, then inspected by eye three separate segments of 
the spectra, away from the NaiD lines. We graded can- 
didates in how well their spectral features matched the 
features of HIP 47513. Four spectra were considered ex- 
cellent matches, and those are the four we chose to in- 
clude in this analysis. 

2.4. X-ray: ROSAT 

We examine X-ray data in order to determine the 
provenance of the SXRB via the X-ray sha dow of the 
LLCC (see Q. The Rontgensatellit (ROSAT; [Truemper 
1992 ) produced a flood of high-quality, all-sky data taken 
during 6 months in 1990 and 1991. Among these data 
were the soft X-ray background (SXRB) maps. We use 
the second data release of these maps, refined from the 
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TABLE 1 

M-DWARFS FROM THE CALIFORNIA PLANET SEARCH 
INCLUDED IN THIS STUDY. 



Star 


V 


B — V 


Date 


Exp. (s) 


GL 250 B 


10.05 


1.42 


2006 Jan 


500 


HIP 47513 


10.38 


1.49 


2006 Dec 


600 


HD 97101 B 


9.95 


1.57 


2005 Feb 


500 


HIP 70865 


10.68 


1.47 


2006 Jan 


500 


HIP 115562 


10.05 


1.47 


2006 Jan 


600 



original reduction to have hig her resolution (12') and fi 



delity (Snowden et al. 1997). The SXRB maps cover 



98% of the sky in 3 wavebands: 1/4 keV, 3/4 keV, and 
1.5 keV. We focus only on the 1/4 keV (C-band) SXRB 
maps in this work, as Hi can act as very effective ab- 
sorber of photons at this energy level, even at the inter- 
mediate column densities of the LLCC. The SXRB maps 
were reduced so as to exclude variable local emission and 
the contribution of point sources. The brightness of the 
SXRB maps are reported in units of 10 -6 counts s _1 
arcmin -2 (Snowdens). 



3. ANALYSIS 
3.1. Radio: GALFA-HI 
3.1.1. HI Profile Saturation and Optical Depth 
Heiles & Troland (20031 measured Hi absorption spec- 



tra and derived Hi excitation temperatures T x for the 
strong radio continuum sources 3C225a, 3C225b, and 
3C237, and found T x = 21.53 ± 1.30, 17.44± 1.80, 13.61 ± 
0.26 K, respectively. These three low temperatures sug- 
gest that the large parts of the LLCC are very cold, < 20 
K, throughout. Given that we find brightness tempera- 
tures in this range for the center of the cloud, we expect 
large parts of the cloud to be optically thick at line cen- 
ter, i.e., saturated. The derivation of accurate Hi column 
densities requires including this line saturation effect; if 
we neglect the saturation, then the derived values are too 
small. 

It is easy to correct for line saturation at the position 
of a sufficiently strong continuum source, because then 
we can measure not only the emission profile but also the 
optical depth profile. However, it is difficult to correct for 
saturation throughout the cloud's area, because without 
a background continuum source our only indication of 
line saturation is the line shape. 

Knapp & Verschuur (1972) attacked this problem by 
assuming that the intrinsic line shape is accurately rep- 
resented by a single saturated Gaussian component and 
used least-squares fits to determine its degree of satura- 
tion and spin temperature T x . However, we have found 
that the stronger Hi profiles tend to be double-peaked, 
so this model is invalid. 

Figure [T] shows an example spectrum towards the core 
of the cloud. Visually, the lines look double. This is 
confirmed by least-squares fitting of three models. The 
solid line is the observed profile, the dotted line the single 
unsaturated Gaussian fit, the dashed line the saturated 
Gaussian fit, and the dash-dot line the double unsatu- 
rated Gaussian fit. The dispersions of the data from 
each fit are indicated on the upper left of each plot. The 
double unsaturated Gaussian fit is the best fit, and we 
find this to be generally true across the highest column 
parts of the cloud. We must therefore consider multi- 
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Fig. 1. — A comparison of a single-component optically thin (dot- 
ted), optically thick (dashed), and a double-component optically 
thin (dash-dotted) fit of the Hi line profile towards a typical, high- 
column part of the LLCC (data are shown by the solid line). Resid- 
uals between the fit and data are overplotted with scale at right. It 
is immediately clear that an unsaturated Gaussian fit does a very 
poor job at fitting the line shape, and that even a saturated Gaus- 
sian has trouble reproducing the double-peaked profile in the data. 
While the double-component optically thin fit largely succeeds, it 
ignores the effects of Hi optical depth, and thus is unphysical. 

component models to accurately fit these data. 

3.1.2. Five Fitting Models 

The best fit in Figure [T] the double unsaturated Gaus- 
sian, neglects the effects of Hi optical depth, and will 
therefore not yield physically reasonable quantities. To 
explore the optical depth of the cloud and its behavior 
over the entire cloud area, we fit five models to each 
pixel's Hi profile: 

(1) : A single unsaturated Gaussian. This is appropriate 

for pixels with weak, single-component lines. This 
model has three unknown parameters: the height, 
center, and width of the Gaussian. 

(2) : A double unsaturated Gaussian. This is appropriate 

for pixels with weak, double-component lines. This 
model has six unknown parameters: the heights, 
centers, and widths of the two Gaussians. 

(3) : A double Gaussian in which the spin temperatures 

are equal. We consider this our most "physi- 
cally reasonable" model; Hi temperatures are de- 
termined by microphysical heating and cooling pro- 
cesses that are not expected to change rapidly with 
position, and thus we do not expect the two com- 
ponents to have different spin temperatures. This 
model has seven parameters: the spin temperature 
and the central optical depths, centers, and widths 
of the two Gaussians. 

(4) : A double Gaussian in which each component has its 

own spin temperature, central optical depth, cen- 
ter, and width, for a total of eight parameters be- 
tween them. An additional (ninth) parameter is 
the line-of-sight ordering of the two components, 
i.e., which lies closer to the observer (and thus ab- 
sorbs the radiation from its more distant sibling). 

(5) : A double Gaussian in which the more distant is op- 

tically thin, and the nearer component is optically 
thick (and thus absorbs the radiation from its more 
distant sibling). This model has seven parameters. 
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3.1.3. Procedure for Fitting Gaussians 

Fitting Gaussians is a difficult and non-unique process 
because the functions are nonlinear with respect to the 
fitted parameters. Such fits require guesses; the equa- 
tions are linearized in a multidimensional Taylor series 
around the guesses and then solved to find corrections 
to the guesses. The corrections are applied and the pro- 
cess is iteratively repeated until convergence or failure. 
The final results depend on the guesses, which makes the 
whole process subjective. See the Appendix for a discus- 
sion of the guessing procedure. 

The fit for model (1) always converges. The fits for 
the other models don't always converge. If a fit does 
converge, we record the mean error of the fitted points, 
which we denote a, and also the mean error of points 
within the velocity ranges where the fitted intensity ex- 
ceeds 0.1 of the peak fitted intensity, which we denote 
aiine- Fits for model (3) usually don't converge because 
for the more distant component the only handle on its 
optical depth is its line shape. We find this is too weak 
a handle, particularly in the presence of noise and the 
line shape modification produced by the nearer compo- 
nent. For this reason, we don't include model (3) in our 
subsequent discussion. 

3.1.4. Model 4 or Model 5? 

Our goal is to settle on three models that apply to the 
whole cloud. Two of these are the single- and double- 
component optically thin models, which are appropriate 
for weak lines, where the effects of optical depth have 
an undetectable effect on our observed line shape. These 
are models (1) and (2). The third, for strong lines, is 
either model (4) or model (5). We choose between these 
by examining their au ne statistics. 

Here we compare au ne values for various model pairs. 
For most of the fits model (2) has a lower au ne than 
model (1) indicating that the vast majority of sightlines 
have double components. 

By doing the same comparison between model (4) ver- 
sus model (5), we find the vast majority of sightlines are 
better fit by model (5). Philosophically, we would have 
much preferred the physically motivated model (4), but 
we must be driven by the data. A second philosophical 
objection is illusory: why should the near side be opti- 
cally thick and the far be optically thin? The answer is 
that the far side may well be optically thick; we cannot 
derive its optical depth because, as discussed in the con- 
text of model (3), all we have is its line shape, and this 
is an insufficient "handle." 

We also make the same comparison between model (5) 
and model (2), i.e., including optical depth effects versus 
assuming optically thin components. We find that while 
a majority of points have lower au ne for model (5), there 
remain a large minority of sight lines that clearly prefer 
model (2) . This is to say that optically thin fits are useful 
along with the optically thick fits. 

For each pixel, we choose the converged model with 
the lowest <Ju ne . Figure [2] shows an image of the cold 
cloud in which the brightness indicates column density 

5 There are schemes to explore the multidimensional parameter 
space for the guesses, but they are computationally intensive and 
not perfectly reliable. 



NOT) [tra ! ] 
10" 



5 21 




ra n 



Fig. 2. — Brightness is <x logioN(HI), and color indicates the 
model chosen for each pixel: model 1 is represented in blue, model 
2 in green, and model 5 in red. The highest column regions are 
typically best fit with model 5, and the lowest column regions are 
best fit with model 1, with intermediate column densities best fit 
by model 2. The position of Regulus is marked by a white star. 

and the color shows the model selected. For each pixel, 
the column density is calculated from the selected model. 

3.2. Infrared 

To generate a map of the IR flux associated with the 
LLCC, first we must remove the bulk of the IR flux which 
is associated with the background Hi, beyond the local 
cavity. We make a map of the background Hi column 
by integrating over the entire Galactic Hi line in the 
GALFA-Hi data and subtra cting the Hi flux associated 
with the LLCC, as found in { 3.1 model (1). Then we fit 



a first-order polynomial to the relationship between the 
background Hi and the IR flux, over the region of the 
map, using a standard least-squares approach. We have 
implicitly assumed that the background Hi is largely op- 
tically thin, typically a good assumption for high latitude 
WNM where T A < T x . We then subtract the expected 
IR background flux from the observed IR flux to deter- 
mine the residual flux we expect to be associated with 
the LLCC itself. The result of this process is shown in 
Figure [3j 

4. CLOUD PROPERTIES 

4.1. Distance 

4.1.1. Lower Limit 

Ruling out Nai absorption along the line of sight to 
HIP 47513 requires detailed analysis for two reasons. 
First, HIP47513's crowded M-dwarf spectrum hides in- 
terstellar absorption features far more eff ectively than 
smoot h A-star spectra like HD 84722 in |Meyer et al. 
(2006). Second, and more problematic, is that HIP 47513 



residual IRAS 100 micron flux 




150 148 
RA [deq] 

Fig. 3. — The residual 100 £im flux from the IRIS maps toward 
the LLCC, shown on a range from -0.7 to 0.7 MJy sr" 1 . The 
LLCC lO 1 ^ cm -2 Hi contour is overlaid in white to guide the eye. 
100 fim flux is quite clearly vi sibl e associated with the LLCC; the 
correlation is shown in Figure 1101 



has a n LSR radial velocity of +4 km s _1 (Nidever et al. 
2002 1 , indistinguishable from the cloud's line-ol-sight ve- 
locity. That is to say, the cold cloud's Nal absorption 
features would fall right the middle of HIP47513's deep 
Nai D lines. 

Because of the difficulty untangling the Na I features of 
the cloud and the star, we compare the Nai D line core 
of HIP 47513 with the line cores of four additional stars 
wit h sim ilar spectra, selected by the method described 
in §2.3[ and shown in Table [I] Our goal is to place an 
upper limit on the Naibetween the Sun and HIP 47513. 

With four stars in hand, we set about normalizing their 
spectra to HIP 47513. In the CPS setup, the Nai D lines 
fall near the edge of one of the HIRES echelle orders. The 
steep blaze function in that region required a somewhat 
sophisticated normalization procedure. Using segments 
of the spectrum 5 A from the Nai line cores, we fit a 
fifth-degree polynomial to both of the Nai line wings. 

With the Na I features normalized across all five spec- 
tra, we zoomed in on the line cores. By eye, it appeared 
that no cold-cloud absorption is present in the HIP 47513 
spectrum (see Figure ffl). We next put a quantitative 
limit on the amount ofNai that could be along the line 
of sight to HIP 47513. 

First, we averaged the four comparison spectra using a 
straightforward mean where each spectrum was weighted 
equally. Then, we subtracted the comparison spectrum 
from the HIP 47513 spectrum and looked for Nai absorp- 
tion features in the residuals. By fitting a Gaussian to 
the residuals at both line cores, we find a maximum ab- 
sorption line intensity of 0.00 ±0.15, equivalent to a Nai 
colu mn density of 0.00 ± 2.1 x 10 10 cm 



In |Meyer et a l. (2006), Nai column is found in 23 stars 
behind the LLCC. We rind that all 8 stars with measured 
cold cloud N m > 2 x 10 19 cur 2 have 7V Na i > 5 x 10 11 
cm~ 2 . As the Hi column towards HIP 47513 is ~ 10 20 
cm~ 2 , the measured 3<r limit of ANai < 6.3 x 10 10 cm~ 2 
implies conclusively that HIP 47513 lies in front of the 
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Fig. 4.— The Nai Di and Nai D 2 line in HIP 47513 and the 
four best-matched CPS stars. The gray lines represent the best- 
matched stars, the black represents their average, and the orange 
represents HIP 47513. The fact that the HIP 47513 lines nearly 
exactly match the average of the lines in the four comparison stars 
indicates there is negligible absorption in Nai Di and Nai D2 from 
the LLCC toward HIP 47513. 



LLCC. 



4.1.2. 



Upper Limit 

As illustrated in Figure [5] our high signal-to-noise 
(S/N « 900) spectrum of Regulus reveals a weak inter- 
stellar Ca 11 K line at an LSR velocity of 3.1 km s _1 . In 
order to give some context for the Ca 11 K line, we have 
also included spectra of the interstellar Nai Di and Ca 
II K absorption toward HD 84722 in Figure [5] HP 84722 
is pos itioned at a distance of 84.8 ± 3.8 pc ( [van Leeuwen 
20071 behind the core of the LLCC. The Na 1 Di spec 



trum of HD 84722 is from Meyer et al.| ( 2006| ) and the Ca 
11 K spectrum is the product of data obtained in March 
2007 with the KPNO coude feed. The latter spectrum 
was reduced in a similar manner as the Regulus spec- 
trum and it represents the sum of six T2KB CCD spectra 
spanning a total exposure time of 11 hours at a velocity 
resolution of 1.5 km s _1 . The measured LSR velocities 
of the LLCC Ca n K and Nai Di absorption toward HD 
84722 are +3.9 and +3.7 km s _1 , respectively. The fact 
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Fig. 6. — An Hi position-velocity slice at <5 = +11°58'30", includ- 
ing the position of Regulus, marked with a white line. The shad- 
ing corresponds to brightness temperature of the Hi, with white 
as higher brightness temperature. The expected velocity (+3.19 
km s _1 ) of the Regulus cloudlet if it were a part of the LLCC is 
marked with a black horizontal line, which is also marked in the 
Regulus spectrum in Figure [5] The position and velocity of the 
Regulus cloudlet is marked with dashed box. The coincidence of 
the Hi velocity of the cloudlet, the expected velocity of the cloudlet, 
and the Ca II K line in Regulus indicate that the LLCC is in front 
of Regulus. 



Fig. 5.— The spectrum of Regulus and HD 84722 at the Ca II 
K line. The spectrum of HD 84722 is also shown for the Na I D 1 
line, as a comparison. The velocity of the Regulus cloudlet, 3.19 
km s~ 1 LSR, is marked with a gray line in the Regulus Ca II K 
spectrum; the near-exact alignment is a clear indication that the Ca 
II K line absorption is created by the Regulus cloudlet. Note that 
the Ca II K line appears broader in the Regulus spectrum because 
of the somewhat lower spectral resolution of the observation. 



that the Nai absorption is much stronger than that of 
Ca II toward HD 84722 is consistent with the character 
of the cold, dense atomic gas at the core of the LLCC. 

Detection of interstellar Ca II absorption towards Reg- 
ulus indicates that some cool interstellar matter must ex- 
ist between us and Regulus. Unfortunately, Regulus is 
not along the line of sight directly toward the body of 
LLCC, but rather about 3° away. Regulus is on the side 
of the LLCC that seems to have a significant amount of 
connected cloud debris (see Figure (21), and indeed Reg- 
ulus is very near a small cloudlet, which we dub the 
Regulus cloudlet. We expect the Regulus cloudlet is a 
component of the LLCC for a number of reasons. First, 
it is only a degree or so away from a long plume of ma- 
terial extending towards higher right ascension from the 
break between the two main LLCC clouds, and therefore 
is most likely a fragment of this plume. Secondly, it has a 
linewidth entirely consistent with the rest of the LLCC. 
Perhaps most tellingly, when we perform a least-squares 
fit to the velocities of the main body of the LLCC with a 
model that accounts for its velocity in 3-space, the veloc- 
ity expected for the position of the Regulus cloudlet were 
it to be moving as a solid body with the LLCC is 3.19 
km s _1 LSR. This velocity is in very nice accord with 
what is observed for the cloudlet (see Figure IBl . This 



also implies that the Regulus cloudlet is largely dynami- 
cally connected to rest of the LLCC, making a scenario in 
which the cloudlet is at a dramatically different distance 
rather contrived. 

The next question is whether the absorption in Ca II 
seen in Regulus is associated with Regulus cloudlet. Reg- 
ulus is on the very margin of the cloudlet, with no ob- 
vious cloudlet CNM Hi detectable directly toward the 
star. The CNM non-detection is consistent with the non- 
detection of Nai typically associated with CNM, as in 
HD 84722 (Figure |5}. Ca II is frequently detected to- 
wards warmer neutral gas or even largely ionized gas, if 
the gas is not so hot as to deplete Ca II into higher ioniza- 
tion states. Therefore, it is not surprising to detect Ca II 
on the edge of an ablating piece of cloud where material 
could easily be heating up and mixing with the warmer, 
more ionized ISM. Perhaps even more conclusively, the 
Ca II absorption in Regulus is at almost exactly the same 
velocity as the nearest part of the Regulus cloudlet (see 



6J. We conclude that we are indeed seeing 
iegulus by the cloudlet, and that therefore 



Figures [5] and 
absorption of 
the LLCC is closer than 24.3 ± 0.2 pc. 

4.2. Temperature 



Crovisier & Kazes (19801 andlHeiles & Troland (20031 



found temperatures for the LLCC ranging from 13 K 
to 22 K for positions towards strong radio continuum 
sources. To determine temperatures for the r est of the 
cloud, we turn to our saturated line fits from 93.11 



4.2.1. Deriving the Spin Temperature for Model 5 

We use the standard radiative transfer equation to de- 
rive the spin temperature T x for the near-side absorb- 
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ing component. Model 5 consists of four independent 
components that contribute to the measured brightness 
temperature Tb ( V) : 

1. The continuum background, which consists of the 
cosmic microwave background (2.8 K) plus the 
Galactic continuum backg round. We e s timat e the 



latter from the 408 MHz |Haslam et al.| ([1982]) sur 
vey by applying a brightness-temperature spectral 
index of —2.6 to their measured value of ~ 15 
K; this gives 0.6 K at 1420 MHz. Thus, we take 
the total continuum brightness temperature to be 



T, 



bkgnd 



3.4 K. 



The contribution from background Hi, which we 
call the WNM contribution because the profiles are 
much wider than those of the cold cloud. This is 
velocity dependent. We represent it by the symbol 
Twnm{Y)- Numerically, it is r epre sented by the 
fitted polynomial mentioned in j]3.1| 



3. The LLCC's more distant optically thin compo- 
nent, which we represent with subscript 0. Its 
brightness temperature in the absence of absorp- 
tion is 



Tb,o(V) = ^oexp 



(V ~ Vo) 2 
2SVJ 



(1) 




Fig. 7. — A histogram of the derived spin te mpe rature for the 
front optically thick component in Model 5 (see j ]4.2| l. The cloud is 
very cold, with modeled t emperatures cons istent with H i absorp- 
tion line observations from |Crovisier &; Kazes| Jl980l and |Heiles fc] 
|Troland| j2003t . 
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4. The closer optically thick component, which we 
represent with subscript 1. Its brightness temper- 
ature in the absence of absorption is 



T B , 1 (V)=T a!>1 (l-exp(-Ti(V))) 



where 



ti(V) = n.o exp 



2 5V? 



(2) 



(3) 
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Fig. 8. — The modeled temperature of the LLCC; color indicates 
the T x i, gray indicates that T Xl i is indeterminate. 



Putting all these together, the emergent brightness 
temperature Tb(V) is 



T B (V)=T BA (V) + 

[Tbkgnd + TwNm(Y) 



(4) 



■T Bfi (V)} exp[-n(V)} 



Our measured profiles arc frequency-switched, which 
means that we measure the quantity Tb(V) — T^kgnd- 
For each pixel, this is the fitted function; we derive the 
spin temperature T x i from the fitted parameters. 

Figure [7] shows the histogram of derived spin temper- 
atures T Xi i for the near component. It peaks near 20 
K and a secondary peak near 26 K, and is well-confined 
within the range 15-30 K. Some anomalously low points 
with very low spin temperatures are presumably the re- 
sult of spurious fits. 

Figure [8] shows the image of T X) i. Regions of high 
column density tend to be warmer. For weak profiles 
with low column density, T x i is often indeterminate, as 
shown by the gray. 



4.3. Distance- Dependent Quantities 

Given the lowe r and uppe r limits on LLCC distance 
derived in E 4.1.1 and [4.1.2 we can determine a range 
for a variety of other quantities that pertain to the cloud. 
We determine the mass by simply integrating the Hi col- 
umn over the area of the cloud, and physical size can 
be determined from angular size. If we make the very 
rough assumption that the cloud is as thick as it is wide 
we can determine a density for the cloud, using a fiducial 
peak column density of 2. 5 x 10 20 cm" 2 . Using the tem- 
perature range found in |4.2[ we can determine a range 
of pressures. We can determine a rough lifetime for the 
cloud by dividing the fiducial scale, i n this c ase the cloud 
width, by the collapse velocity (see { 4.3.1). Results are 
shown in Table 2. 

Our assumption that the cloud's thickness along the 
line of sight is comparable to its width (i.e. the cloud 
is tubul a r), is at odds with the result from 
Troland (2003) that t his cloud is thin along t 



Heiles fc 
le line of 

sight (i.e. sheet- like). |Heiles &: Troland (2003) assumed 
the cloud was outside the local cavity, and thus 4 to 9 
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Fig. 9. — A map of the velocity difference between the two col- 
liding components. Brightness is tx (TAV) ' 5 , summed over both 
velocity components. Color indicates the velocity difference be- 
tween the two components. 

times more distant than we now know it to be. This dis- 
tance assumption gave an overestimate of the width of 
the cloud and thus an o verestimate of its aspect ratio. 
Heiles & Troland (2003) also measured the thickness at 
low-column (N H ~ 3 x 10 19 cm' 2 ) sight-lines to 3C225a 
and 3C225b, about 10 times lower column than the thick- 
est part of the cloud in the present analysis. This low col- 
umn density gave a lower expected thickness for a given 
measured temperature and assumed pressure, further ex- 
acerbating the inconsistency with our work. If the cloud 
were 10 times flatter along the line of sight than we have 
assumed, the derived pressures would be 10 times higher 



far out of t he standard ISM pressure range (e.g., Wolfire 



et al. 



the c 



2003 ) . This result validates our present picture ot 



oud as more tubular than sheet-like. 



4.3.1. Kinematics of the Two Velocity Components 

We found that two velocity components usually pro- 
vide better Hi line fits than a single velocity component. 
Model (5) allows us to discriminate between the veloci- 
ties of the two components, because the optically thick 
foreground component is closer. Accordingly, comparing 
their velocities allows us to determine whether the com- 
ponents are approaching or receding from each other. 

There is an overwhelming tendency for the front com- 
ponent to have higher velocities than the rear. This indi- 
cates that the components are approaching one another, 
typically with a velocity difference ~ 0.4 km s _1 . 

Figure [9] images the velocity difference between the 
two components, with color indicating velocity difference 
and lightness the total profile area, i.e., the sum of the 
two component areas. We only show pixels which had a 
successful two-component fit, where the two components 
were not degenerate. Generally, the edges of the cloud 
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Fig. 10. — The relationship between Hi column and infrared flux 
for the LL CC u nder an Hi optically thin assumption (contours) and 
the fits in {j3~l"1 (grayscale). The dust-to-gas ratio of 48 X 10 -22 MJy 
sr — 1 cm 2 for the j ]3.l| fits is shown as a black and white dashed 
line, while the standard high-latitude dust-to-gas ratio found by 
|Schlegel et al.| | |1998) of 67 xlO -22 MJy sr _1 cm 2 is shown in a 
black dotted line. INote that the optically thin assumption leads to 
a strange dependence of the dust-to-gas on the Hi column. 

(where the intensity is low and the line-of-sight thickness 
is small) have larger relative velocities. In these portions 
of the cloud, the cloud is getting thinner with time by 
~ 0.4 pc Myr -1 . If we assume the two components are 
not more separated than the width of the cloud, we find 
the collision time to be approximately 1 Myr. 



TABLE 2 

Quantities for the LLCC and LRCC based on the 

DETERMINED DISTANCE RANGE. DENSITY, PRESSURE, AND 
COLLISION TIME MEASUREMENTS DEPEND ON THE ASSUMPTION 
THAT THE CLOUD IS AS THICK AS IT IS WIDE. 



quantity 



range 



unit 



distance 
LLCC Hi mass 
LLCC length 
LLCC width 
LRCC length 
density 
pressure 
collision time 



11.3 - 24.3 
0.235 - 1.07 

2.8 - 5.9 
0.25 - 0.54 
13 - 28 

320 - 150 
9600 - 2250 

0.6 - 1.3 



pc 

pc 
pc 
pc 
cm -3 
K cm" 3 
Myr 



4.4. Dust Content 

To determine the large-grain dust content of the cloud 
we investigate the relationship between the Hi column 
density and the 100 /urn emission. Figure [TU] shows the 
relationship between these two quantities over the cloud 
area. It is clear from this plot that there is a very nar- 
row and tight linear relationship between Hi column and 
IR flux. This serves as an important confi rmation that 
the fitting procedures used in Section |3.1| to determine 
the saturated Hi column densities are not wildly incor- 
rect. For comparison we show the relationship between 
the IR flux and the Hi colum n u nder an optically thin 
assumption, model (1) (Figure 10 contours). These two 
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quantities are clearly not linearly correlated, confirming 
that an optically thin assumption for the Hi in the LLCC 
is wrong. 

The cloud has a 100 ^m to Hi column value of 48 
xlO -22 MJ y sr - 1 cm 2 . This is somewhat lower than the 



result from Schlegel et al. (1998) for low column density 
areas of 67 xlO -2 ^ MJy sr _i cm 2 . It is also lower than 
the standard Galactic values found in Boulangcr fe Per-| 
ault| ( |l988[ ), about 100 xl0~ 22 MJy sr" 1 cm z , although 
consistent with their measurement towards Auriga. Sim- 
ilarly, it is lower than t he average value f ound for Galac- 
tic clouds measured in Peek et al. (2009), but consistent 
with the values found for the LI, L7, and L8 clouds found 
in that work. We conclude that this cloud has either a 
lower-than-expected overall dust grain density, or a pop- 
ulation that has somewhat larger grains than typical for 
the ISM that maintain a somewhat lower temperature. 
In either case the cloud is only marginally anomalous in 
IR production via grains. 

5. IMPLICATIONS FOR THE LOCAL CAVITY 

Authors since Verschuur's discovery paper have re- 
marked upon the observational exceptionality of the 
LLCC: it is remarkably cold, bright, and long in the 
Galactic Hi sky. The cloud is also exceptional spatially; 
it is one of just a few clouds with column densities above 
10"" 



20 



et al. 



in 2 within a few dozen pc of the sun (Vergely 
20101 RL08). This makes it a excellent and rare 



test element for studying the contents and dynamics of 
the gas that surrounds us within the local cavity. 

5.1. Interaction with the CLIC 

In recent years much pioneering work has been done 
on the closest known interstellar clouds to the sun (e.g 
Redfield fc Linsky||2008l |Frisch|[2008l |Redfield fe Falcon| 
2008p . The core results from these investigations is that 
the sun resides within the CLIC, a collection of clouds 
within about 15 pc of the sun, and that these CLIC 
clouds move as solid bodies past the sun in roughly the 
same direction. These clouds have typical column densi- 
ties near 10 1 
authors 



cm , with temperature near 10 4 K. Many 



Audit fc H enncbcllc 2005} |Vazquez-Semadem 
BT Heitsch et al. 2006 ) "nave conjectured that 
Ids are formed by the convergence of flows of 
RL08 pointed out that the LLCC is 



et al. 20i 
CM M clou, 
warm neutral gas 
coincident on the sky with the faster moving Gem cloud 
and the slower LIC, Auriga, and Leo CLIC clouds, and 
suggest that the collision of these clouds may be respon- 
sible for the formation of the LLCC. 

Our new lower limit on the LLCC of 11.26 ± 0.21 pc 
seems to rule this out. All of the clouds suggested by 
RL08 to be in collision to form the LLCC are detected to 
be, at least in part, closer than 11.1 pc (the LIC envelops 
the sun; Gem is detected within 6.7 pc, Leo within 11.1 
pc, and Aur within 3.5 pc). The CLIC clouds certainly 
span some range of distances, and the extent to which 
they are entangled is still largely unknown; it is possible 
that the LLCC exists at some more complex interface 
between these clouds. That said, the simple picture in 
which the LLCC is at the interface of two independent 
colliding clouds seems unrealistic given the new distance 
lower limit on the LLCC. 

While the distance of the LLCC does seem to separate 
it somewhat from the CLIC, the velocity of the LLCC 




Fig. 11. — A diagram showing the heliocentric velocities of the 
15 clouds described in RL08 (circles) and the velocity of the LRCC 
described in HaudlO (cross). The velocities associated with the 4 
clouds closest on the sky to the LRCC are outlined in a dashed gray 
line. The LRCC has a very similar velocity to the RL08 clouds, 
suggesting they are associated to each other or some other parent 
cloud or population. 



may imply that the CLIC and the LLCC are more re- 
lated. HaudlO fit the LRCC (the ribbon of clouds that 
contains the LLCC) with a traversing, rotating, expand- 
ing, ring. It is unclear if this model has physical signif- 
icance, as such structures aren't typically found in dy- 
namic models of the ISM, but the results are striking. 
The overall velocity of the ring is 37.4 km s _1 toward 
I = 213°, b = -11°. 1 in the HaudlO fit (when converted 
to the heliocentric reference frame), which is nicely con- 
sistent with the velocities of the CLIC clouds found in 
RL08 (Figure 11). It is possible that this consistency is 



a coincidence, out it seems more likely that the CLIC 
clouds and the LRCC are related to each other in some 
way, or to a parent cloud population. 

5.2. Shadowing of Soft X-rays by the LLCC 

As we outlined in the Introduction, there is significant 
controversy over the origin of soft X-rays seen coming 
from all directions in the sky. In the local hot bubble 
theory, soft X-rays at low and intermediate latitudes em- 
anate from a pervasive hot gas filling the local cavity, 
with the excess of soft X-rays at high latitude coming in 
part from the halo. In the SWCX theory, the isotropic 
component of soft X-rays com es from the very nearby 
(~ 100 AU; " ' 



505} 
SM. 



Stone et al. 20051) interaction of the solar 



wind with the surrounding ISM. We can differentiate be- 
tween these theories by looking for absorption of these 
X-rays by the LLCC. 
In the direction of the LLCC the wall of the local cav- 



ity is between 100 pc and 150 pc away (Meyer et al. 
20061, and we have demonstrated th at t he LLCC is be- 
tween 11.26 pc and 24.3 pc away (J 4.1). We have de- 
termined the Hi column of the LLCC rather carefully 
([3.1 ), and the local bubble wall has a column density 
of ~ 3 x 10 20 cm~ 2 in the direction to the LLCC. For a 
reasonable spectral energy distribution over the 1 /4 keV 
band of ROSAT, an optical depth of 1 corresponds to an 



Hi column of 1.05 x 10 2U cm" 2 (Snowden et al. 1997) 
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Fig. 12.— The RASS 1/4 keV (C-band) brightness as a function 
of LLCC transparency in the C-band, given the measured Hi col- 
umn. Diamonds represent pixels in the RASS map, and the line 
represents the fit found in Equation IB] The red bars represent 1 a 
errors in Ic, foreground a n d Ic, total- The fact that 398 ±38 Snow- 
dens emanate from in front of the cloud rules out a n model with 
no SWCX and a smoothly emitting LHB (see i ]5.2.1| l. 



To determine the amount of flux coming from in front of 
the LLCC we fit the equation 

In (a S) = In t . + /^, , , p -T C (HI(a,S)) 

J C \ LX ^ U ) — 1 C ,/ oreground i ^Cbackground^- 

(5) 

where Ic (a, 5) is the flux detected toward the cloud in 
the SXRB 1/4 keV (C-band) map, and t c (HI (a, 6)) is 
the optical depth to C-band X-rays caused by the LLCC 
as a function of position on the cloud. Ic, foreground 
and Ic, background are the unknown X-ray flux from in 
front of and behind the cloud respectively. To do this 
we use a standard least-squares fit. The standard er- 
rors derived from a least-squares fit are correct under the 
assumption that the deviation from the model is Gaus- 
sian and uncorrelated. Neither of these assumptions is 
true in our case, since there are obvious structures in 
the SXRB 1/4 keV maps in this direction that are far 
from random, and will generate strong, unmodelcd er- 
rors. To account for this we use the disp l aceme nt map- 
ping technique described in Peek et al. fj2009 ), which 
determines the errors by examining the fluctuation in 
the result as the cloud is placed in different positions 
on the nearby sky. We then use this fluctuation to 
better estimate the errors. Using this method we find 

Ic, total = Ic, background + Ic, foreground = 708 ±11 SnOW- 

dens and Ic foreground = 398 ± 38 Snowdens (see Figure 
12b. 



To constrain models for the origins of the soft X-rays, 
we combine the above result with the results toward 
the interme diate latitude (b = —34.0 4) molecular cloud 



MBM 12 in Kuntz et al.| ( |1997| ) and |Kuntz fc Snowden 
( 2000[ ). In this analysis the^uthors found that MBM 
2 was at the back of the local bubble, as it shows very 
little, if an y, shadowing in t he C-band, is at a distance 
of ~ 90 pc (|Sfeir et al.|1999|), and with a C-band flux of 



347_ 10 . Any model we construct for the origin of the C- 
band X-rays must account for these two measurements. 

5.2.1. Model A: Local Hot Bubble and Halo Emission 



In this model the LHB is filled with a homogeneous 
gas such that the intensity, Ic, is simply proportional to 
the path length of LHB gas. We add to this some un- 
absorbed halo emission from the halo , equivalent to the 
trans- absorptive emission discussed in |Kuntz fc Snowden| 
(2000). We can show easily that model A is insufficient 
to explain the data. We find that 389 ± 38 Snowdens of 
soft X-ray flux are not attenuated at all by the LLCC, 
and thus are generated closer to us than the LLCC, along 
a path length no longer than 25 pc. If these X-rays were 
generated by a hot, local cavity-filling gas, we would ex- 
pect at most 708 Snowdens *(25 pc / 100 pc) = 177 
Snowdens to originate from in front of the cloud, as the 
ratio of path lengths is no more than 25 pc / 100 pc. 
Adding any halo emission only exacerbates the incon- 
gruence. This expected value of 177 Snowdens is incon- 
sistent with our result by more than 5 cr, thus strictly 
ruling out model A. 

5.2.2. Model B: Solar Wind Charge Exchange and Halo 
Emission 

In model B we assume the LHB does not contribute 
at all to the soft X-ray background. Instead we appeal 
to solar wind charge exchange, which we model as a sin- 
gle isotropic, all-sky inten sity in the C-band. We note 
that models of SWCX by (|Kout roumpa et al.l|2008| ) find 
an expected maximum dipolar variation of ~25% in the 
all-sky C-band maps from the direction of motion of the 
heliopause through the LIC, although they do not model 
the outer heliospheric region, which may largely mitigate 
this variation. An all-sky value of 347 Snowdens for the 
SWCX is 1.1 a off from our measurement of the C-band 
intensity in front of the LLCC, and also satisfies the in- 
tensity measurement toward MBM 12. The background 
intensity off the LLCC of 708 - 389 = 319 Snowdens 
is explained by the halo emission. We note that if the 
background flux co uld as easily be described by the "Hot 
Top" framework of |Welsh fc Shelton| ( |2009[ ) as by halo 
emission. In either case, model B is consistent with these 
results. 

5.2.3. Model C: Local Hot Bubble, Solar Wind Charge 
Exchange, and Halo Emission 

In model C we allow all three components to contribute 
to the soft X-ray background. Since model B fits the data 
successfully, we concern ourselves with determining the 
maximal LHB contribution that is consistent with the 
data. To have a non-zero contribution from the LHB, 
the true value towards the LLCC must be less than the 
true value toward MBM 12 (as MBM 12 is more distant), 
which is 1.1 cr from our measurement. If we take 3 a as 
a credible limit, we find that the LHB has a maximal 
emission of 1.1 Snowdens / pc, abo ut 4 times lower than 
found in |Kuntz fc Snowden (2000), with 247 Snowdens 
coming from SWCX. At 2 <r, the emission is only 0.52 
Snowdens / pc, with 300 Snowdens coming from SWCX. 

6. CONCLUSION 

We have examined the the LLCC using Hi emission 
data, diffuse far-infrared dust emission data, optical stel- 
lar absorption line data, and soft X-ray background data. 
We have concluded that the LLCC is within 24.3 pc (Reg- 
ulus), but beyond 11.3 pc (HIP 47513). Through careful 
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fitting of the Hi line, we have found that the cloud is typ- 
ically composed of two components, moving toward each 
other at 0.4 km s" 1 . We also have found typical tempera- 
tures for the LLCC between 15 K and 30 K, with column 
densities ranging up to 2.5 x 10 20 cm" 2 . We have found 
that the LLCC is too far away to be created by the colli- 
sions of the CLIC clouds discussed in RL08, but that the 
velocity determined for the parent LRCC cloud is very 
consistent with the CLIC clouds, pointing to a common 
origin. Perhaps most importantly, we have shown that 
the extent to which the LLCC absorbs C-band radiation 
is inconsistent with a homogeneous Local Hot Bubble in- 
terpretation for the isotropic soft X-ray background. We 
have shown that this absorption instead favors the solar 
wind charge exchange interpretation, though we do not 
strictly rule out a model in which both contribute the 
the soft X-ray background. 

While we have answered many questions in this work, 
we have also posed many question. Is the multi-peaked 
profile a standard aspect of ultra-cold CNM? Is there 
significant variability in the SXRB emission within 25 
pc? What is the chemical composition of the LRCC? We 
hope to address these questions with new observations 



of the LRCC. We intend to complete the map of the 
LRCC above 6 = -1° with the GALFA-Hi survey, using 
archival and scheduled observations. These observations 
will allow us to further constrain the SXRB emanating 
from within the distance to the cloud, as some regions 
appear have high column density similar to the LLCC. 
The larger cloud area to examine may also allow us to 
narrow the distance range with new and archival high 
resolution stellar absorption observations. Observations 
are planned of various metal lines in the optical and UV 
towards stars behind the LRCC, which will allow us to 
better understand the physical state of the cloud, includ- 
ing its metallicity, ionization fraction, and pressure. 
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APPENDIX 

We handled the "guessing problem" mentioned in §3.1.3| using the following sequence for each pixel: 

1. We first fit model (1), the unsaturated single Gaussian. We considered the velocity range -4.0 to +10.0 km 
s" 1 and defined the "line" velocity range as -1.2 to +7.5 km s . To define the guesses, we first subtracted a 
third-order polynomial fit to the off-line profile values. We took the height guess equal to the maximum profile 
value within the line velocity range. We took the guesses for center, and width (our widths are always FWHMs) 
equal to the first and second moments of the HI profile in the line range. Finally, using the results of this fit, 
we fit simultaneously for the polynomial and the Gaussian parameters over the full velocity range; this fit was 
almost always successful, and when it wasn't we used the results from the first fit. 

2. We use the solutions from step 1 to generate guesses for model (2). Guesses for the two centers are equal to the 
model (1) fitted center ±X, where X = 0.3 times the model (1) fitted width; for the widths, we always used 0.8 
km s _1 ; and both heights were equal to model (1) fitted height. 

3. We use these solutions from step 2 to generate guesses for model (3). The guesses for centers and widths are 
equal to the solutions, the guesses for optical depths are both equal to 1.5, and the guesses for spin temperature 
are adjusted to make the unabsorbed peak brightness temperature of the components' emission lines equal to 
their model (2) peak brightness temperatures. The guesses for model (5) are obtained in a similar fashion. 

4. The guesses for models (4) are obtained in a similar fashion as for model (3), except that for T x and r we 
performed fits over a closely spaced grid of T x , so no guessing was involved for these two parameters. 
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